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We argue that the anomalous power asymmetry observed in the cosmic microwave background 
(CMB) may have originated in a cosmic bounce preceding inflation. In loop quantum cosmology 
(LQC) the big bang singularity is generically replaced by a bounce due to quantum gravitational 
effects. We compute the spectrum of inflationary non-Gaussianity and show that strong correlation 
between observable scales and modes with longer (super-horizon) wavelength arise as a consequence 
of the evolution of perturbations across the LQC bounce. These correlations are strongly scale 
dependent and induce a dipole-dominated modulation on large angular scales in the CMB, in agree¬ 
ment with observations. 
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I. INTRODUCTION 

The Planck team has reported evidence of anoma¬ 
lous features in the large-scale cosmic microwave back¬ 
ground (CMB) which point towards violations of statis¬ 
tical isotropy [T]. Some of these anomalies were already 
observed by WMAP, and have been now re-confirmed 
[2], hence reducing the possibility of instrumental origin 
or systematics. Although the associated statistical sig¬ 
nificances are still inconclusive, the intriguing possibil¬ 
ity that these features are traces of fundamental physics 
beyond the simplest inflationary models have attracted 
considerable attention in the theoretical community. 

Among the different anomalies detected, the power 
asymmetry —unequal power in different regions of the 
sky—has received significant attention, and it is also the 
main target of this paper. This asymmetry was first mod¬ 
eled by adding a modulating factor to an otherwise sta¬ 
tistically isotropic temperature distribution [3] 

6T{h) = (l + A(n)) 5Tiso{n) (1) 

where Aiji) — J2 lm ^LMYLAiifi) is the modulating 
function. Planck has reported no evidence for Alm dif¬ 
ferent from zero for L > 2, but a statistically significant 
(about 3 (t) L = 1, dipolar modulation has been found 
[T]. Furthermore, the effect of this dipole on the angu¬ 
lar power spectrum Ci has been detected only for low 
multipoles with i < 64. More precisely, after separat¬ 
ing the £-range into bins of size = 64, only the first 
bin shows a non-vanishing signal, with average amplitude 
Ap^i = 0.07 ± 0.02. Observations are compatible with 
Ai=i = 0 for larger values of £. This scale dependence 
implies that the simple parametrization 0 is insufficient 
to account for the observed modulation. The theoretical 
challenge is therefore to find a mechanism able to produce 
a dipolar modulation present only on large scales, with 
negligible contribution to the quadrupole, octupole, etc, 
and furthermore respecting the existing constraints on 
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the CMB: a remarkably Gaussian, almost-scale invariant 
spectrum of adiabatic perturbations. Not surprisingly, it 
has been difficult to find a completely satisfactory model 

Hi- 

One of the first and most compelling ideas to gener¬ 
ate such a power asymmetry was introduced by Erick- 
cek, Carroll, and Kamionkowski in It relies on the 
fact that the presence of a very long wavelength, super- 
Hubble perturbation of a curvaton field will induce a 
dipole modulation in the observed spectrum, provided its 
wavelength is still short enough that its amplitude shows 
a gradient across the sky. This long mode could originate, 
for instance, as a remnant of the pre-inflationary epoch. 
This idea has been implemented by several authors and 
improved in different directions (see e.g. [SHE]). 

Another interesting proposal, introduced in [9], works 
instead in a single-field inflationary model, and generates 
the power modulation through a non-Gaussian coupling 
between observable scales in the CMB and even larger, 
super-Hubble scales. The reason why non-Gaussianity 
can induce anisotropies in the observed spectrum, even 
if the underlying statistics is isotropic, is simply because 
a typical realization looks significantly more anisotropic 
if the underlying distribution is non-Gaussian. Or in 
other words, because observable modes couple to the par¬ 
ticular realization of the long wavelength modes in our 
Universe, which is generically anisotropic. This model, 
therefore, requires a mechanism to generate large cor¬ 
relations between very different scales. The consistency 
relation proposed in m tell us that it is difficult to find 
a realistic single field model producing that type of non- 
Gaussianity. 

In this work we introduce a scenario which has features 
in common with both of the previous ideas. We consider 
a single field inflationary model preceded by a bounce 
described by loop quantum cosmology (LQC). This sce¬ 
nario has been analyzed in great detail, both for the back¬ 
ground space-time and for perturbations (see pTVH4] for 
reviews). In short, the evolution of perturbations across 
the bounce excites quanta out of an initial vacuum, and 
as a consequence the onset of inflation is reached in an 
excited state, rather than in the Bunch-Davies vacuum. 
The presence of those perturbations, remnants of the 
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pre-inflationary phase, have an important impact on the 
non-Gaussianity generated during inflation. We compute 
those non-Gaussianity and the modulation they produce 
on observable scales, and show that the observed power 
asymmetry can have an origin on the quantum bounce 
preceding inflation, while still respecting all observational 
constraints. We set c=l but keep G and h explicitly. Nu¬ 
merical values are given in Planck units. 


II. LQC AND THE POWER SPECTRUM 


In LQC, the mean effective space-time geometry is de¬ 
scribed by equations which incorporate the leading quan¬ 
tum corrections to general relativity. For instance, the 
Friedmann equation reads [TTl [T^ 
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where H is the Hubble rate, p the energy density, and 
Pmax ~ 0.4/5pi its Upper bound, which is a fraction of the 
Planck energy density ppi. In this paper the matter sec¬ 
tor is assumed to be a single scalar field with mass m. 
More complicated potentials can be incorporated. How¬ 
ever, the effects we are seeking have a quantum grav¬ 
ity origin, and the form inflaton potential produces sub¬ 
leading contributions. At low energies p <C Pmax, Eq. 
([^ becomes indistinguishable from the general relativis¬ 
tic Friedmann equation. However, close to the Planck 
scale quantum gravitational effects break the linear re¬ 
lation between and p, and bring H to zero while p 
attaints its maximum value. This is a quantum bounce 
that ‘bridges’ a contracting and an expanding phase of 
the universe. 

The theory of cosmological perturbations in LQG that 
we use in this paper was developed in [16], and has 
been used in [igiiHi to study the evolution of pertur¬ 
bations across the bounce through the end of inflation, 
and to compute the primordial power spectrum of tensor 
and scalar perturbations. This scenario is based on first 
principles, and trans-Planckian issues can be addressed 
squarely. The free parameters relevant for phenomeno¬ 
logical studies are the value of the inflaton field at some 
arbitrary reference time, e.g. the bounce, 4){tB) '■= 4 >b, 
its mass m, and the initial state of perturbations at some 
initial time. Refs. [laiisi have explored the predictions 
for the power spectrum across the parameter space and 
contrasted the results with observations. 

From the viewpoint of inflation, the effect of the pre¬ 
ceding bounce translates to an excited state |/3) for per¬ 
turbations at the onset of slow-roll. Therefore, the LQC 
pre-inflationary evolution can be conveniently encoded in 
the Bogoliubov coefficients and /3fc relating |/3) and the 
Bunch-Davies vacuum, i.e. relating the mode functions 
of curvature perturbations TZk{t) that result from the 
pre-inflationary evolution, and the Bunch-Davies modes 
72.®^ (t). These coefficients can computed numerically for 
a choice of the free parameters. For (pp = 1.22 and 


m = 1.10 X 10“®, both in Planck units, and “Minkowski- 
like” vacuum initial condition^ for perturbations at one 
Planck second before the bounce, we obtain that the av¬ 
erage number of quanta present at the onset of inflation, 
given by |/3fep, is approximately 10“^ for the reference 
mode fc* which today corresponds to 0.002Mpc“^. For 
the longest wavelength mode we can directly measure, 
fcmin « fc*/8.9, we obtain |/3fc„i„P = 1.2, and |/3fcp ~ 1/k 
for lower values of k [HI- Other choices of initial data 
give similar results (see [nidH] for details). 


The inflationary scalar power spectrum in the presence 
of an excited state |/3) is given by Pn{k) = P^{k) \ak + 
/IfeP, where P^{k) is the Bunch-Davies result. Fig. 
shows the numerically computed Pnik) for the choice m 
parameters mentioned above. Two energy scales play an 
important role in the power spectrum. First, LQC intro¬ 
duces a new energy scale kpQc/a{tB) '■= ^/RbJQ « 3.21, 
that is directly related to the space-time scalar curvature 
at the bounce Rb = dSTrpmax ~ 62. A second scale is 
provided by the value of curvature at the onset of acceler¬ 
ated expansion kj/a{ti) := \fRjJ% ~ 10“®, which occurs 
at time t/—notice that ti is not the onset of slow-roll, 
which happens at later times. Modes k > kpQc are ‘in¬ 
side’ the curvature radius—i.e. its wavelength is shorter 
than the curvature radius—during the bounce, and exit 
during the slow-roll era. These perturbations reach the 
onset of slow-roll in the Bunch-Davies vacuum, and their 
power spectrum is indistinguishable from standard re¬ 
sults. Modes ki < k < kpQC exit the curvature radius 
soon before the bounce and re-enter after it, to exit again 
during slow-roll. The crossing process around the bounce 
time amplifies the amplitude of those modes, which then 
reach inflation in an excited states. Finally, modes k < kj 
are outside the Hubble radius during the entire evolution, 
even at the onset of inflation, and their spectrum is sig¬ 
nificantly suppressed. The reference mode /c*, that today 
corresponds to 0.002Mpc“^, is approximately one third 
of kpQc, fc* = 0-86 k^Qc- Gonsequently, the LQG cor¬ 
rections to the observable power spectrum only appear 
for the largest angular scales in the sky, that correspond 
to multiples ^ 30, and are only significant for the low¬ 
est multiples. Furthermore, the LQG corrections signif¬ 
icantly increases the power for modes with wavelengths 
larger than the Hubble radius today. This power en¬ 
hancement conveniently reaches a maximum around kj, 
at which perturbation theory is still well under control. 


^ As explained in m, although Minkowski-like initial conditions 
do not produce a state with the desired ultraviolet behavior— 
i.e. a Hadamard or fourth adiabatic order state—one can always 
modify the initial data for sufficiently ultraviolet modes to make 
the state ultraviolet-regular. Since such modification do not af¬ 
fect observable predictions, we do not describe the details in this 
paper, which can be found in dsiiia 
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IV. THE BISPECTRUM 

At leading order in perturbation theory, non- 
Gaussianity is characterized by the bispectrum 
B']i(ki,k 2 ,k'i), defined in terms of the three-point 
function: 


“ (27r)^(5^(fci + k2 + ^ 3 ) B-ji{ki, ^ 3 ) ■ 


FIG. 1: Scalar power spectrum at the end of inflation for (j)B = 
1.22 and m = 1.10 x 10“®, and “Minkowski-like” vacuum 
initial conditions at one Planck second before the bounce [18] . 
Gray points show the numerical result for individual modes. 
The spectrum is rapidly oscillatory, and its average is shown in 
black. The spectrum is amplified for low wave-numbers ki < 
k < kLQC, but the enhancement is only significant for the 
very low multipoles in the GMB, for which the observational 
error bars are large because cosmic variance. Therefore, this 
power spectrum is compatible with current observations. 


III. BACK-REACTION ON THE 
INFLATIONARY GEOMETRY 


An important question is whether the energy in the 
quanta present at the onset of slow-roll produces signifi¬ 
cant back-reaction on the inflationary space-time geom¬ 
etry which cannot be neglected. The energy density on 
the perturbations can be computed and compared to the 
vacuum result. The difference in energy density between 
the state |/3) and the Bunch-Davies vacuum can be writ¬ 
ten as 


(3) 


where the integral in k is extended from —oo to oo. It 
is well known that the difference of expectation values of 
the energy-momentum tensor between states that are at 
least of fourth adiabatic order is always finite, hence no 
renormalization is required in (|^. Numerical evaluation 
shows that the ratio of Ap with the background energy 
density po, is Ap/po « 10“^ at the onset of slow-roll, 
and decreases exponentially with cosmic time. Therefore, 
since the energy density in the Bunch-Davies vacuum is 
known to be negligible small, the back-reaction of the 
state 1/3) on the inflationary geometry can be neglected. 


The non-Gaussianity generated from excited states in in¬ 
flation has been analyzed by several authors (see e.g. 
[T^22j b It was pointed out in m that the main charac¬ 
teristic of the associated bispectrum is an enhancement 
in ‘squeezed’ configurations which involve very different 
scales, fci « ^2 S> ^ 3 . These non-Gaussian correlations 
originate from quantum interactions between particles 
that are present at the onset of slow-roll. In realistic 
models, including the scenario presented in this paper, 
the conclusions of the consistency relation proposed in 
m are still satisfied in the limit ki/k^ —)■ 0, but never¬ 
theless there are important effects for small, but finite, 
values of ki/k^. 

The full expression for the bispectrum B']i{ki,k 2 ,kff) 
as a function of the state |/3) was obtained in 
and is given in the appendix. We have used those re¬ 
sults to numerically computed B'fi{ki,k 2 ,k 3 ) at the end 
of inflation. The main difference with previous analysis 
is that the state |/3) is now computed from a quantum 
gravity framework, rather than postulated. The result¬ 
ing bispectrum, therefore, carries information about the 
pre-infiationary evolution. Those effects are more clearly 
displayed by plotting the ratio with the standard Bunch- 
Davies bispectrum B^{ki,k 2 ,kffj (see Fig. [^. As ex¬ 
pected, the Bispectrum is strongly peaked in squeezed 
triangles, where the ratio B-ji/B^ shows a growth pro- 

portional to —i.e. proportional to product of the 

number of quanta in the modes fci and k^, as expected for 
an enhancement originated from particle interactions. As 
for the power spectrum, the enhancement only appears 
for modes kj < k < kLQc- Therefore, there are strong 
correlations between the largest wavelengths we can di¬ 
rectly observe (k « k^qc) and super-Hubble modes with 
k values in the range kj < k < kLQc- The bispectrum 
becomes negligibly small when the three momenta are 
in observable scales, hence respecting current constraints 
on non-Gaussianity. This scale dependence of the bispec¬ 
trum will play an important role in the next section. 

Recall that the three modes involved in bispectrum 
must form a triangle, so it is really a function of two 
momenta, e.g. B'fi{ki, ks). The shape of the bispectrum, 
shown in Fig. can be understood by writing its domi¬ 
nant contribution in the squeezed limit x := k^/ki <C 1: 
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B-ji K. Ac IS.Tiiki) X Re ft 

h 


1 _ pikirio 1 _ pik^rio 

/i7r^^+/2 


X _ gikt-no 
1 + (1 + /i) X 

1 _ gifc3'?0 


(1 + m): 


(1-Xi); 


1 + (—1 + /i) X- 


(4) 


where = ki ■ ks, kt = ki + k 2 + k^, ki = kt — 2fci, and 
A-jzik) = (^) evaluated at Hubble-crossing 

during inflation. The functions ft contain the informa¬ 
tion of the state |/3) and are defined in the appendix. 
The parameter ryo is the value of the conformal time at 
the onset of slow-roll. For the example (pB = 1-22 and 
m = 1.1 X 10“® we obtain ryo^* ~ 10^, so all observable 
modes are deeply inside the Hubble radius at the onset 
of slow-roll. 

The largest contribution to Q comes from very 
squeezed configurations for which k^rjQ <C 1. In that limit 
B-fi is dominated by the terms proportional to ft and /s. 
This contrasts with the example considered in [20l [22] of 
a scale invariant excited state |/3), and where additionally 
the analysis was restricted to observable scales for which 
Tjoki 3> 1 for all i. In that case the leading contributions 
come from the terms proportional to fi and / 2 . 



FIG. 2: Ratio of the inflationary bispectrum for the excited 
initial state arising from the LQC pre-inflationary evolution 
versus the Bunch-Davies bispectrum, as a function of k 2 and 
fcs for a fixed ki — 0.22fc*. The plot shows the range of 
k 2 and ks allowed by the triangle condition ki + k 2 + ks = 
0. The bispectrum B-jz is highly scale dependent and shows 
a prominent enhancement for squeezed configurations fcs ^ 
k2 ~ ki. 


It is important to note that the non-Gaussianity com¬ 
puted here are generated during slow-roll inflation. Ex¬ 
tra contributions to the bispectrum will certainly arise 
from the evolution across the bounce. However, since 
the state of perturbations is very close to vacuum at the 
bounce time, those contributions are expected to be sub¬ 
dominant. This expectation is indeed borne out in ex¬ 
plicit computations in bouncing models p5l IM] , where it 
is shown that the non-Gaussianity generated across the 
bounce are much smaller than the bispectrum shown in 
Fig.H particularly for squeezed triangles, which are the 
relevant configurations for this paper. It is therefore rea¬ 


sonable to assume that the leading contributions to the 
bispectrum in squeezed configurations are dominated by 
inflationary non-Gaussianity. 


V. NON-GAUSSIAN MODULATION 


We follow some of the ideas presented in Ref. [3] to 
compute the CMB power modulation arising from cou¬ 
pling with super-Hubble scales. The statistics of the clas¬ 
sical Bardeen potential $ for observable scales is modified 
in the presence of a given long wavelength perturbation 
A?{ki) if the spectrum is non-Gaussian. Concretely, the 
two-point function acquires off-diagonal contributions of 
the form [53 j^Sj 


=P$(fc) i27Tfd{k + k') + Gik,ki)<A?{ki) (5) 


where ki = —{k + k') (i.e. it closes a triangle with k 
and fc'), P<s>{k) is the <I>-power spectrum, and we take 
G{k,ki) = ^BTiik,ki)[An{k)Aniki)]-^. For bispectra 
peaked in squeezed configurations the leading contribu¬ 
tions to G{k,ki) come from the regime ki k, for which 
k K, k'. Therefore, the correlation function is close to be 
diagonal. The off-diagonal terms break both homogene¬ 
ity and isotropy. These terms vanish if we average over 
as it must be, since our model respects these sym¬ 
metries at the fundamental level. But in the particular 
realization of <I>^^ chosen by our Universe, the modula¬ 
tion term may be important, and may produce deviations 
from homogeneity and isotropy in the CMB much larger 
than would be expected from a typical realization of a 
Gaussian spectrum. 

The off-diagonal terms in (HI source analogous terms in 
the covariance matrix of temperature spherical harmonic 
coefficients 

{aimOfim') = ^uAmra'Gl + ^ A^M Q-mrn'M + Cb) , 

LM 


where is a Wigner 3-j symbol, and Gt is the 

standard angular power spectrum. The momentum de¬ 
pendence of the kernel G{k,ki) factorizes for the domi¬ 
nant contributions in a power law for k (except for typical 
oscillations which are unimportant for our computations 
and can be averaged). The kernel can be expanded in 
Legendre Polynomials 

G(fc,fc;) = PLifk), 


where fi = k ■ ki, and gL{ki) encodes the ki dependence. 
Then, the mean value of the modulation amplitude, av¬ 
eraged over M, for a typical realization of the long wave¬ 
length modes can be computed as (see [S] for further 
details of the computation) 


{Am 


Gg{—aL) 

2GeiO) 


dki kf 
(27r)3 


1/2 


\9L{ki)\^Pm 


( 6 ) 
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FIG. 3: Amplitude {Al{£)) as a function of £ for L=l,2,3, for parameter values (pB = 1-22 and m = 1.10 x 10 and 
“Minkowski-like” vacuum initial conditions at one Planck second before the bounce. 


where the coefficients Ci,{aL) have been defined as the 
CMB temperature power spectrum computed by replac¬ 
ing ris —>■ Us + Q!l- Therefore, if aL = 0 for all L, i.e. if 
G(fc, ki) does not depend on k, the modulation amplitude 
is ^independent. 

The value of {Al{£)) for different multipoles L can 
be computed numerically from the bispectrum shown in 
Fig. and some of its properties can be qualitative 
understood already from equation The leading con¬ 
tribution comes from the terms proportional to ft and 
/a, which are dominated by a dipole L = 1. The terms 
proportional to /i and /2 contribute to a subdominant 
quadrupole. The octopole L = 3 is subdominant com¬ 
pared to the quadrupole, and this hierarchy continues 
for higher multipoles. Since the factors ft in Q scale 
approximately as {kik)~^ for every z, we have « 1 
for all L, and therefore the amplitude A{£) is expected 
to decrease with £. 

Figurej^shows {Al{£)) for L = 1, 2,3 as a function of £. 
The amplitudes {Al{£)) are all scale dependent, and such 
that (Ai) > 10 {A 2 ) > 10 (A 3 ). The average value of (Ai) 
for € < 64 is in agreement with the observed value = 
0.07 ± 0.02. Other values of our free parameters (fs and 
m, and of the initial conditions for perturbations, can 
decrease or increase the value of (Ai). For instance, we 
find that choosing vacuum initial data for perturbations 
far into the past of the bounce, significantly decreases the 
amplitudes of all multipoles. 

There has been some debate based on symmetry ar¬ 
guments about whether a statistically homogeneous and 
isotropic bispectrum can generate a dipole modulation. 
Our computation is an example where the answer is in the 
affirmative. Furthermore, as shown in m. a careful anal¬ 
ysis reveals that symmetry arguments restrict A^+i to be 
suppressed with respect to A^ for even L. The hierarchy 
of multipoles we find here, (Ag) > (Ai) > (A 2 ) > ...., is 
therefore in agreement with those restrictions. 

VI. DISCUSSION 

We have presented an inflationary scenario in which 
perturbations start the slow-roll phase in an excited 
state, rather than the Bunch-Davies vacuum. This state 


arises from the pre-inflationary evolution provided by 
loop quantum cosmology, in which the big bang singu¬ 
larity is replaced by a bounce. Two new scales appear 
in the problem, kj and k^qc, related to the onset of the 
exponential expansion and the bounce, respectively. The 
number density of quanta at the onset of slow-roll is sig¬ 
nificant only for the range ki "^k k^qc- During infla¬ 
tion these excitations induce non-Gaussian correlations 
which we have computed. To the best of our knowledge, 
this is the first computation of non-Gaussianity in LQC. 
The result is compatible with existing observational con¬ 
straints. Furthermore, large correlations arise between 
the longest modes we can observe, with k « krqc^ and 
super-Hubble modes with k > kj. We have shown that 
those non-Gaussian correlations, which involve super- 
Hubble modes, are able to modify the observed power 
spectrum at large scales, inducing correlations between 
CMB angular multipoles ^ 30 that differ in A£ « 1. 
These correlations are strongly scale dependent, and pro¬ 
duce a power asymmetry in agreement with observations. 

Other observed anomalies at large scales [T] —parity vi¬ 
olation, power suppression for a bunch multipoles around 
I = 20 , multipole alignments, etc.—seem also to require 
correlations between low multipoles qualitatively similar 
to the ones obtained in this paper. Whether the agree¬ 
ment is also quantitative will be analyzed in future work. 

We emphasize that the fact that the power spec¬ 
trum shown in Fig. [^is amplified for low wave-numbers 
ki < k < kLqc, is not necessarily in conflict with obser¬ 
vations since: i) the enhancement is only significant for 
the very low multipoles in the CMB, for which the obser¬ 
vational error bars are large because of cosmic variance; 
ii) the observed power suppression is only significant for 
a few multipoles around £ = 20 , indicating that the ef¬ 
fect is more likely to originate from correlations between 
multipoles, rather than a suppression of the primordial 
two-point function. 

Some of the features appearing in our scenario— 
remnants, large amplitude perturbations associated to 
an infrared scale kj, correlations with super-Hubble 
modes, etc.—have been identified in previous phe¬ 
nomenological analysis (see e.g. O [9]) as ingredients 
needed to account for some of the observed anomalies at 
large scales. Here, these features arise from a concrete 
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quantum gravity proposal based on first principles. 
Therefore, our results provide further motivation to 
consider the observed anomalies as real physical fea¬ 
tures, which have origin beyond the simplest inflationary 
models, rather than statistical flukes or instrumental 
noise. Future work will be addressed to provide further 
robustness to the model introduced here and to extend 
its quantitative predictions. 
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Appendix A: Inflationary bispectrum from an 
excited state 


The expression for the scalar bispectrum generated 
during slow-roll inflation when the state of perturbations 
1/3) is given by a Bogoliubov transformation of the Bunch- 
Davies vacuum, with coefficient ak and /3fe, is given by 

[HHH] 


B 


n 


(/ci, ^ 2 , fca) — 2 — 2ry -I- e ^ ^ ^ -I- 


(Al) 


4e 


kl 


Re 


1 _ piktrjo 1 _ pikirto 1 _ pik^-qo l _ pik^qo 

ft - , - +/l-^- +/2-~- +/3-^- 

kl k2 k^ 


I -I- 2 cyclic permut. fci —?> ^2 —>■ fca , 


where kt = ki + k 2 + k 3 , ki = kt — 2ki, and the parameter 
r]o is the value of the conformal time at the onset of slow- 
roll. The functions fi contain the information of the state 

m-- 


ft=[l + F{ki) (1 -b F{k 2 )) + cyclic perm.] , 


and /2 and /a can be obtained by cyclicly permuting 
the momenta in /i, and F{k) = |/3fcp -b a^Pk- The 
star indicates complex conjugation. Also, — 

(Ay evaluated at Hubble exit during infla¬ 

tion. 


/i = F(fci)*[l + F{k2) + F(fc3)] - F{k2yF{k ^), 
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